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We discuss the status and some ongoing upgrades of the ZFITTER program for applications at e^e~ colliders 
LEP 1/SLC, LEP 2, GigaZ, and TESLA. The inclusion of top quark pair production is under work. 



1. Introduction 

One of the great achievements in theoreti- 
cal and experimental physics over the last three 
decades is the unique calculation of quantum cor- 
rections in the Standard Model (SM) of fundamen- 
tal particle interactions, following the pioneering 
work of G. 't Hooft and M. Vehman [ |l|, |, |, 
and their observational evidence at the 10 a con- 
fidence level [ ^ . Especially the electroweak sec- 
tor of the SM is being investigated at the high 
energy e^e~ colliding facilities LEP 1 and SLC, 
and now at LEP 2, with unprecedented precision. 
This, though, could well be surpassed in the fu- 
ture e.g. by the GigaZ mode of the TESLA Linear 
ColHder (LC) project [|, §. One main focus, be- 
sides the search for direct signals of New Physics, 
was and still is on studying the properties of the 
neutral and charged weak gauge bosons. 

Here we concentrate on the measurement of 
different cross sections and asymmetries in fer- 
mion pair production. Fermion pair production 
is an essential tool to determine parameters like 
the mass and total and partial decay widths of 
the weak neutral gauge boson Z and the neu- 
tral current fermionic couplings. The theoreti- 
cal break-through for this was the first complete 
evaluation of electroweak quantum corrections to 
e+e" M^M^ in the Weinberg model in [ 
with a collection of the scalar one-loop integrals 
first given in [ || . 

For today's predictions at the few per mil level. 



^Talk given at "Loops and Legs 2000" , Bastei, Germany, 
April 9-14, to appear in the Proceedings. 



this calculation can only be seen as the first es- 
sential step. It did not include such important 
ingredients as a refined treatment of the Z boson 
resonance, of QCD and higher order corrections, 
or realistic QED hard bremsstrahlung effects. 

During the last twenty years many dedicated 
efforts were invested to perform the necessary im- 
provements, to implement them into computer 
codes which may be applied in the analysis of 
experimental data and in regularly repeated com- 
parisons of the codes. A representative collection 
of the underlying expressions, their implementa- 
tions and numerical comparisons may be found in 
[ |, |l| ; see also [|T|, ^ for individ- 

ual comparisons. The most recent workshop on 
LEP 2 physics is yet under way [ |l7| . Prominent 
examples of such codes are: 

- ALIBABA, 

- BHM, 

- KORALZ, KK, 

- TOPAZO, 

- ZFITTER. 

On the latter, the semi-analytical Fortran pro- 
gram ZFITTER [ we wiU now focus. EarHer 
program descriptions are [ |2^. The core of 
ZFITTER relies on a complete electroweak one- 
loop calculation [ |2^, It is also 
based on many additional formulae; from [ 18 
one may extract a list of papers that ZFITTER 



uses in addition, [ 26, 27, 2g, 29, 30, 31, 32, 33 
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2. Status of two-fermion codes at LEP and 
higher energies 

ZFITTER [|l^ calculates radiative corrections to 
the muon decay constanU |2|, H, ||, |ll[, to tlie 
Z decay width [ ||, ||, |l|l , and to the W decay 
width [ ^1^. Cross sections and asymmetries are 
treated in a semi-analytical approach. Improved 
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containing the complete virtual weak and QCD 
corrections are convoluted with analytical flux 
functions p{s' / s) (radiators) for the QED correc- 
tions [ 1^ |2|, H, 111. This is done in a one- 
dimensional numerical integration over the final- 
state invariant mass squared s' = M'^j [ [l8| . 
The especial Bhabha scattering case is treated 
in an effective Born approximation only follow- 
ing [ |Tl|, |3^. As an example, the initial-state 
corrections to total cross sections for s-channel 
processes are reproduced in short (v = \ ~ s' / s): 
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The radiators p™*"'"-'^™, and p™* with S and 
Ht are determined in [ ^ , and including all rel- 
evant higher order terms e.g. in [ |l8j. Eq. (|l|) 
can be straightforwardly generalized to different 
asymmetries Aps, ^poi, ^l_r etc. or to scattering 
angle distributions like da/dcos9, then with dif- 
ferent effective Born terms and radiators. Kine- 
matical cuts to the final-state phase space may 
also be applied [ |8|, l23l El ~ 



p5| , |95[ . Also a more 
model-independent description of cross section 
observables, e.g. in form of an S-matrix approach, 
can be thought of [ H, ||, |9|, || ||, |8||44| . 

Concerning the status of present codes for 
LEP 1 and SLC applications, one can summarize 

§ A Fortran bug was corrected in ZFITTER v. 6. 30 [ ^ 
resulting in a 0.3% shift of the W partial widths. 



that the level of precision for cross section predic- 
tions in s-channel fermion pair production is now 
better than 10^** on the Z boson resonance and 
better than 0.3 x 10""^ for center-of-mass energies 
y/s = Mz ± 3 GeV (single contributions). For 
the recently updated branch with cuts on final- 
state maximum acoUinearity and minimum ener- 
gies in ZFITTER [ |6[ |5), this is demonstrated in 
Table |l|, comparing with program TOPAZO [ Q 
when the QED initial-final state interference is 
switched on. For earlier comparisons with acoUi- 
nearity cuts also consult [ p7| p8|, p9[ . 

At LEP 2 center-of-mass energies up to roughly 
200 GeV an agreement of the codes at the or- 
der of few per mil is being achieved [ |l^ . This 
was studied earlier for the two-fermion programs 
ZFITTER [ ^ |l8| BHM [|lo|, 0, TOPAZO [|§, and 



102| for diffe rent cuts on min- 
104 



KORALZ/KK 

imum s' [ |ll], |103| , |14| , |15| , |1U4|| or on maximum 
acoUinearity and minimum energies of the final- 
state fermions [ |6[ Similar studies 
of B habha s catte ring include the code ALIBABA [ 
|105|| in [ ^ 106 . The conclusion was that pre- 



cisions are of the order of few per mil to 1% for 
LEP 1 or LEP 2 energies respectively, excluding 
a radiative return to the Z by sufficiently strong 
cuts. Compared with the final experimental pre- 
cisions at LEP 1/SLC and LEP 2 these theoreti- 
cal accuracies are satisfactory [ 107 . It is worth 



mentioning that for the envisaged accuracy sev- 
eral recent improvements had to be undertaken, 
e.g. a better treatment of initial- and final-state 
pair production corrections, the exponentiation of 
initial-final state corrections, convolution of the 
ZZ and WW box corrections, etc. As a result, 
ZFITTER v.6.30 [ |l8| became more complicated 
and also slower, and the same might be true also 
for the other codes. 

Thinking of applications of the ZFITTER code 
and of other programs at a future e^e~ Linear 
Collider (LC) with much higher luminosities and 
energies, the above observations at LEP/SLC en- 
ergies might only prove to be preliminary and a 
further upgrade could become necessary. Two dif- 
ferent scenarios may be envisaged for the use of 
the ZFITTER code at a LC: 

First, precision physics could again be per- 
formed on the Z boson resonance, but then with 
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Table 1 

A comparison of predictions from ZFITTER v. 6. 30 [ |l8) ( Jun 2000) and TDPAZO v.4.4 [ Q for muonic cross 
sections and forward-backward asymmetries around the Z peak. First row is without initial-final state 
interference, second row with, third row the relative/ absolute effect of that interference in per mil (Af^ = 
91.1871 ± 0.0021 GeV, = 125 GeV, toj = 173.8 GeV, ag = 0.119, Aa(f2j(^z) = 0.0280398089). 
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0.15909 




+2.75 
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1000 times the luminosity of LEP 1 [ |, |] in order 
to search e.g. for virtual effects from a SM or MSSM 
(minimal supersymmetric) Higgs boson or from 
supersymmetric particles [ 108, ^. With such a 



GigaZ option, experimental accuracies could in- 
crease by a factor of 100 or more [|,|]. Updates 
of the codes could include higher order QED ra- 
diative effects, effects by beamstrahlung, or an 
update of the still critical Bhabha scattering case 
when demanding high precisions. What is more 
complete electroweak two loop calculations [ [109 



might become necessary; see also the many talks 
on this topic at this conference. At energies up t o 
roughly 800 GeV like for the Tesla project [ |no| , 



issues like experimental and theoretical precisions 
are still quite vague, but it is clear that the de- 
mands on codes for two-fermion production will 
be quite challenging due to higher beam energies, 
higher luminosities, and improved analysis tech- 
niques. On the one hand, electroweak and QED 
corrections become equally important which may 
demand a critical look at the numerical validity 
of the usually applied improved Born approxima- 
tion at higher energies. On the other hand, higher 



order electroweak corrections will also grow in im- 
portance with increasing energies [ 111|. First 



studies of codes ZFITTER, TOPAZO, and KK show 
some evidence that an agreement of 5 per mil to 
1 per cent can be reached for the complete energy 
|ll| (see Fig. |). 
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3. Top pair production at LC energies 

At a future LC, top pair production will be 
studied in detail. One of the main topics of in- 
vestigation will be the analysis of the elemen- 
tary properties of the top quark, including the 
measurement of its mass toj, total and partial 
decay widths, and couplings to the weak gauge 
bosons. The recent progress in calculations at 
the tt threshold itself is summarized e.g. in [ 



114 and is continuously being updated in [ 115 



A determination of the weak neutral vector and 
axial-vector couplings of the top quark to the Z 
boson will be performed in the perturbative re- 
gion sufficiently above ^/s « 350 GeV by mea- 
suring cross sections, forward-backward and po- 
larization asymmetries. Here ZFITTER may natu- 



4 



e* e —> /i* /u. {n y) 
0° < i5 < 1 80° 



QED: ISR + FSR, no IFI 




zf: zfitter v.6.22 

TZ: TOPAZO v.4.4 
KK: KK2f v.4.12 



*'^Vs7s>0.9164 



Vs'/s> 0,8397 
Vs'/s>0.6441 



a(TZ)/CT(ZF) 



0.99 



CT(KK)/a(ZF): Vs'/s>0.9164 
□ A 



NO INI PP INI PP EEX: LLO(a')+0(a') CZZX: O(a') 

rii" 



2M« 



2M, 



V s (GeV) 



Figure 1. Cross section ratios for muon-pair pro- 
duction with different s' cuts for codes ZFITTER 
v.6.22 (Oct 1999) [ TOPAZO v.4.4 [ |§, KK 
v.4.12 [ [l0| from 60 to 800 GeV cm. energy; 
without initial-final state interference (INI PP: 
initial-state pair production; LL: leading logarith- 



mic terms) [112 



rally step in to do the job. The one key point now 
to be considered is the inclusion of the mass of the 
final-state heavy fermions. Up to now final-state 
masses could be neglected for applications at the 
LEP or SLC center-of-mass energies. 

For the massive case, the complete one-loop 
electroweak corrections in the SM with the run- 
ning of the QED coupling, fermionic self energies, 
vertex corrections, and electroweak box contribu- 
tions have been determined in [ 116 . A recent up- 



date of the QCD effects was given in [ 117 . Work 



is right now in progress treating the electroweak 
effects in form of effective couplings together with 
form factors in an improved Born approach [see 
term <7^{s') in (|l|)], suitable for a quick calcula- 
tion of interesting observables with the ZFITTER 



code [ 118 . For results within the MSSM please 
refer to [ [119^ pO| . 

In the remaining part of this contribution, we 
focus on the top mass dependent radiative cor- 
rections from QED bremsstrahlung which corre- 
spond to the radiators yO^(s'/s), a = ini, fin, int, 
A = T, FB [see (||)] . One may expect at most 
some additional several per cent photonic correc- 
tions from the finite top quark mass. Initial-state 
corrections to ctt [ 121 1 and final-state effects to 
the angular distribution da /d cos (i.e. also to 
Afb) are known [ |l2|, |3^, pT^ . 

For the massive cross sections ctt,^fb = 
aT / ap-B, da/dcosd, analytic expressions like (|l|) 
would be nice to have. Soft and virtual QED cor- 
rections with final-state masses have been known 



for a lon g time for pure QED [ 123 and also in 
119 1 . For the hard massive radi- 



116 



the SM ^ ^ , ^ 

ators p'\{s' / s; mj / s) the situation is the follow- 
ing: With an s' cut only, hard photonic correc- 
tions for total cross sections ar^s) with mf ^ 
were shown in [ 121 1 and [ |3^ for the initial- 
and final-state corrections respectively. Having 
in mind a cut on the cosine of the scattering an- 
gle cos I? of one final-state fermion, the differen- 
tial cross section da/dcosO first has to be deter- 
mined. For this case, only formulae with final- 
state bremsstrahlung are given in [ |3^; with- 
out cuts, the total cross section and forward- 
backward asymmetry had been presented earlier 
in [|12|. 

It should be mentioned that the QED correc- 
tions to cross section contributions ctl and an for 
a left- or right-handed polarized beam are not 
automatically given by the results for unpolar- 
ized cross sections; there are additional radia tors 
that do depend on the initial helicities [ 124 1. Of 



course, the effective Born results also have to 
be replaced by the corresponding massive terms 
cr° _R with polarization. Nevertheless, the predic- 
tions for A^fi, as being measured at SLC at the Z 
resonance using a polarized e~ beam, may be to a 
very good precision described by the massless ra- 
diator functions (as is presently done in ZFITTER) 
since there hard bremsstrahlung is strongly sup- 
pressed, and, in addition, for polarization asym- 
metries the photonic corrections cancel to a large 
extent (see e.g. [B3]). 
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3.1. The distribution dafp^/ds' 

The hard-photon contribution from QED inter- 
ference to total cross sections for heavy fermions, 



e-(pi) + e+{p2) -> f{p3) + Kpa) + 7(fc), 



(4) 



will be calculated here by two methods; one is 
tensor integration. The final-state phase space 
is split into two two-particle phase spaces. The 
(//) rest frame is boosted with respect to the 
center-of-mass system (cms) . The phase space in- 
tegration is then conveniently carried out in two 
steps: First, the final-state tensor T^'^ is inte- 
grated over p3 and pa in the (//) rest frame, 
leaving Lorentz covariant expressions in k'' and 
qP = p^ + p1. Contracting then with the initial- 
state tensor I^iy, the remaining integration over 
the cosine of the photon polar angle with respect 
to the beam axis is carried out in the center-of- 
mass system. Evidently, the cms fermion produc- 
tion angle is not accessible in this approach and 
thus neither the angular distribution nor ApB- 
The integrated distribution gets: 



ds' 



E 



(5) 



where 



ds 



-{Z'Z) = -Aa' Q,Qf 



(6) 



1 — s'/s 



2m} 



L 



and 

Pf 
Xzis) 



•1+^ 

,1-/3/ 



72 27ra s - M| -h iMzTz ' 



(7) 

(8) 
(9) 



For Z-boson exchange we have: Vf 

9 f 

2Q f sin 9w and af = with /| 



r - 

^3 



= -1/2, 



for example, we have simply — Q f , af = 0, etc. 
and X7 1- Massive contributions are propor- 
tional to 2m'} / s ■ Lf which drop out for (3f ^ 1 
in the massless case. The formula (^) is quite 
compact. 

It may of course also be used in the con- 
text of searches for extra heavy gauge bosons 
Z' e.g. through ZZ' mixing effects (see e.g. [ 
|7|, |12|| and program ZEFIT [Isl). The mass- 



less limit is given in Eq. (1.3.14) in 

Numerical results are shown in FigTg. We ob- 
tained that figure also with another approach (see 
the next section), where we also include a short 
discussion. 



3.2. 



The distribution da^pj^/ds' 



Another phase space parameterization makes 
the fermion production angle accessible. The ap- 
proach goes back to [ ^ and is explained in sec- 
tions (1.2) and (1.5.1) of [ |l8|; see also [ ^ for 
more details, where also the first analytical inte- 
gration is described. After that, the integration 
is much simplified, and it is here where presently 
the Fortran program topfit.f [ 124| starts the 
numerical treatment. We have some understand- 
ing of the resulting integrand, and it seems to 
us that the remaining analytical integrations may 
also be performed while retaining a finite final fer- 
This would open the way to 



mion mass m 



f 



124] 



Qe = —1. For the contributions from 7-exchange, 



treat the massive initial-final state interference on 
the same footing as the other corrections, namely 
with only one numerical integration (over s'). 

The numerical effects from final mass correc- 
tions to the initial-final state interference QED 
bremstrahlung arc depicted in Fig.||. There, both 
the invariant mass distributions da^p^/d^/s' are 
shown without and with final-state mass effects 
for TO/ TOt = 174 GeV and = 500 GeV. 
The tree-level cross section, shown for compari- 
son, is cr° = 0.51 pb. 

For large invariant masses, \/s' ^/s, the in- 
frared peak clearly dominates. It has to be regu- 
larized by soft and virtual corrections. The main 
difference between the massive and massless cases 
is a suppression of cross sections near the thresh- 
old y/s' « 2TOt. This suppression is stronger for 
da'^pg than for dalp* . The total corrections from 
bremsstrahlung are roughly by a factor four larger 
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Figure 2. Invariant mass cross section distribu- 
tions daip^/dVs' and da'j^ydVs'. 



for dapB than for dar in both cases. Including 
final-state masses mt, however, reduces the radia- 
tive corrections for both dapB and dar again by 
roughly a factor of four. 

4. Conclusions 

Summarizing, the semi-analytical Fortran pro- 
gram ZFITTER for e^e~ — s- // with radiative 
corrections is applied in electroweak precision 
physics at LEP and SLC. With its quick cross 
section evaluations and precision of about 10^'^ 
at the Z peak and of order few per mil to 1% up 
to ^/s ~ 800 GeV, the code is well-equipped for 
the latest data analyses at present accelerators 
and to become an important tool for data-fitting 
in fermion pair production at a future LC. 

In the massive case, the code has to be up- 
graded for the still missing branch of continuum 
top quark pair production at a LC. First attempts 
in this direction have been presented here. With 



the upgraded code, or a condensed version topf it 
for continuum tt physics, a determination of the 
weak couplings of the top quark may be stud- 
ied. ZFITTER, or topf it respectively, would then 
describe all fermion pair final-states including ra- 
diative corrections and different realistic cuts for 
the complete range of energies e.g. at DAPHNE, 
LEP/SLC, and at a LC, i.e. for ^/s « 0(1 GeV) 
up to 0(1 TeV). However, at DAPHNE energies 
a specialized code might be better suited in view 
of the different hadronic final states there which 



are observed besides fermion pairs [ 126, 127 1 
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